Two-state quantum walk on two- and three-dimensional lattices 
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We present a new scheme for a discrete-time quantum walk on two- and three-dimensional lattices 
using a two-state particle. We use different Pauli basis as translational eigestates for different axis 
and show that the coin operation, which is necessary for one-dimensional walk is not a necessary 
requirement for two- and three- dimensional walk but can serve as an additional resource. Using this 
scheme, the probability distribution from Grover walk using four-state particle and other equivalent 
schemes on a square lattice using coin operation is reproduced. We also present the Hamiltonian form 
of evolution which can serve as a general framework to simulate, control, and study the dynamics 
in different physical systems. 
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I. INTRODUCTION 



Quantum walks, developed as the quantum analog 
of the classical random walks p~H4] first emerged as a 
powerful tool in the development of quantum algo- 
rithms [6H2] • Subsequently, its rich dynamics is consti- 
tuting as a framework to understand and simulate the 
dynamics in various systems. For example, they have 
been used to explain phenomena such as the break- 
down of an electric-field driven system [10] and mecha- 
nism of wavelike energy transfer within photosynthetic 
systems [TTJ [12] , to demonstrate the coherent quan- 
tum control over atoms [13] and localization of Bose- 
Einstein condensates in optical lattice [14] and to ex- 
plore topological phases [15] . The quantum walk evolu- 
tion are widely studied in two forms: continuous-time 
quantum walk (CTQW) [16] and discrete-time quantum 
walk (DTQW)gJ UZHH]. In this letter we focus on 
the DTQW evolution which is defined on the position 
Hilbert space, H p and the coin (particle) Hilbert space, 
T-L c . During last few years, experimental implementa- 
tion of the DTQW has been demonstrated with en- 
ergy levels in NMR[23, ions[24j[25], photons [261J29] . 
and atoms [30] . These experimental implementations us- 
ing one-dimensional (ID) DTQW model on a two-level 
system using a degree two coin operation have opened 
up a new dimension to simulate quantum dynamics in 
physical systems like the recent demonstration of local- 
ization of photon's wavepacket [3Tj . Now the immedi- 
ate interest would be to extend the implementation to 
two-dimensional (2D) and three-dimensional (3D) lattice 
structure with the available resources. This will give way 
to simulate and explore the possibility of mimicking the 
dynamics in various naturally occurring physical systems. 
One of the extension to the 2D is the Grover walk which 
is defined on a four-level particle with a specific initial 
state [32j [33] . An alternative extension to higher (d) di- 
mensions is to use a d coupled qubits to describe the inter- 
nal states [34j [35] . This is extremely challenging with the 
available resources to implement it experimentally. To 



overcome this challenge, an alternative 2D DTQW us- 
ing a two-state particle was very recently proposed. By 
evolving the particle in superposition of position space 
in one dimension followed by the evolution in the other 
direction using Hadamard coin operation was show to be 
equivalent to four-state Grover walk [36] . 

In this letter we present a new scheme using different 
Pauli basis as translational eigenstate in different axis 
and show that the two-state particle walk can be imple- 
mented on a physically relevant 2D, square, triangular, 
kagome, and 3D, cubic lattice structures. Using basis 
vectors of the three Pauli matrices is very common in 
quantum optics experiments and various other physical 
systems making this scheme implementable in the present 
experimental setups. We also present Hamiltonian form 
of the evolution which can serve as a general framework 
to simulate, control, and study the dynamics in differ- 
ent 2D and 3D physical systems. We then show that 
the coin operation which is required for ID DTQW is 
not a necessary requirement for 2D and 3D DTQW but 
can be used as an additional degree of freedom to con- 
trol the dynamics. In particular, we demonstrate that 
the probability distribution obtained using Grover walk 
and alternative two-state walk on a square lattice [36] 
is effectively reproduced without a coin operation in our 
scheme. This further reduces the resource required for 
the experimental implementation. 



II. ONE-DIMENSIONAL DTQW AND ITS 
HAMILTONIAN FORM 



The standard form of DTQW evolution on a two- 
state particle in ID lattice is defined on a coin (c) and 
the position (p) Hilbert space H = He ® H p . The 
basis states of 1~L C are the internal states of the parti- 



cle, | 1} 



and | t) 



and they are also the 



eigenstates of the Pauli matrix os 



The ba- 
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sis states of % p is described in terms of \ip z ), where 
z G I, the set of integers associated with the lattice sites. 
Each step evolution of ID DTQW is described using a 



2 



quantum coin operation B a3 (0) 



cos(0) sin(0)l, 
-sin(0) cos(0)\ m 
which evolves the particle (coin) into the superposition 
of the basis states followed by the unitary shift operator 

s* 3 = EJU)UI®l^-i)(^l + lt)(tl®l^+i)(^|], 

which shift the state of the particle in superposition of 
the position space. Therefore, the effective operation for 
each step of the DTQW is written in the form: 



<?3 

The state after the 
is given by, |* t ) = 



1}. 



(1) 



t step evolution of the DTQW 
[WasWItfin), where \% n ) = 
(cos(£/2)| \) + e l7] sin(<ty2)| t>) ® |^o), is the initial state 
of the particle at position z = 0. The coin parameter 6 
controls the variance of the probability distribution of the 
walk and 9 ^ 0,7r/2 is required to spread the amplitude 
in superposition of position space. 

The net evolution of each step of the DTQW in the 
form of Eq. ([!]), can also be generated by the time inde- 
pendent effective Hamiltonian H a3 (0). That is, W a3 (0) = 
e -iH a3 (9)r w h ere fi=\ and r is the time required to im- 
plement one step of the walk. To obtain an expression 
for H(6) we will expand Eq. ([!]) and rewrite to the form: 



cos{6)e- iPzlT sm(0)e- iPzlT 
- sm(9)e +iPzlT cos(0)e +iPzlT 



(2) 



where P z is the momentum operator whose action on 
all position space in the Z axis is local such that 
e ±%Pzlr \^z) = \ip z ±i) [U [5], r is the time required to 
implement each step of walk length I (lattice separa- 
tion). Hereafter, we will consider the time required to 



implement the unit length of each step of walk to be 
one (r = I = 1). To obtain the Hamiltonian form, 
Wcr 3 (6) = e~ lH(7 z(6) is written as 



iH a3 (0) = In 



cos(6)e- iPz sm(0)e- iPz 
sm(6)e+ iPz cos(6)e+ iPz 

\n(A) = V\n(X a3 )V-\ 



x z 
.0 A+ 



(3) 



(4) 



is the diagonal matrix composed of the eigenvalues 



= cos(0) cos(P z ) =F a/cos 2 ((9) cos 2 (P z ) - 1 (5) 



of matrix A where \+\ z =\ z \+ = l. 



V = 



and 



v- 1 



sin((9)e if 
1 



sin(#)e* 



sin(<9)e ijF 
1 



A+— cos(#)e* 



(6) 



2v/cos 2 (6>) cos 2 (P z )-l 2 v /cos 2 (6')cos 2 (Pz)-l 
-sin(0)e +iP z -A_+cos(g)e^g 
_2 v /cos 2 (6>) cos 2 (P z )-l 2 v /cos 2 (6>)cos 2 (P z )-l 



(7) 



are the matrix composed of eigenvectors of A and its 
inverse. By substituting these elements into Eq. ([3]), 



In 



2Wcos 2 (<9)cos 2 (P z ) - 1 



cos(#) s'm(P z ) 
sin(6>)[sin(P z ) - icos(P z )} 



- sin(0) [sm(P z ) + i cos(P z ) 
cos(#) sin(P^) 



cr 3 



(8) 



where sin ( P z )\^ z ) 
cos(P z )|^) = \{\^ z 



-1> + I^+1». 



|^+i)) and 



III. TWO-STATE PARTICLE WALK ON 
DIFFERENT LATTICES 

A. Square Lattice and Cubic Lattice 

A simple 2D lattice structure is a square lattice with 
four direction of pro pagation and two quantization axis, 
X and Z [Fig.Fga) 



Similarly, a simple 3D lattice is a 
cubic lattice with six direction of pr opaga tion and three 
quantization axis, X, Y", and Z [Fig. |l(b)"] . 



Two-state particle DTQW on a square and cubic lat- 
tice can be realized by quantizing the evolution using 
different Pauli basis states as translational eigenstate for 
each axis in the lattice structure. For a walk in ID (Z 
axis) we used basis states, | \) = |+) crs and | t) = |— ) a3 
of Pauli operator cr 3 as translational state. Similarly, for 
X and Y axis we will use the basis states, 



\+)*2 



1 

i 



; \-U - 

5 I )cr 2 



1 

i 



1 



1 



(9) 
(10) 



of Pauli operators <j\ 



1 

1 



and <t 2 



-i 

1 



as trans- 
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FIG. 1: (a) Square lattice with two direction of propagation in each 
X and Y axis, (b) Cubic lattice with two direction of propagation 
in each X, Y and Z axis. Evolution in X, Y and Z axis are 
quantized by the basis states of the Pauli operators ai, o<i and (T3, 
respectively. 



lational states. The choice of a particular Pauli basis for 
particular axis is purely conventional and all the matri- 
ces hereafter are represented in the basis formed by the 
eigenvectors of cr 3 . The general form of the coin opera- 
tion in any of the basis states |±) CTc , of the Pauli operators 
<j a with a = 1, 2, 3 will be 

B aa (0) = cos(0)|+) CTq (+| + sin(0)|+> ffa H 

-sm(e)\-) aa (+\ + CO s(e)\-) aa (-|. (11) 

The shift operator for each axis of the square lattice [Z 
and X) and the cubic lattice (Z, X, Y) will be 



,(+| ® \lpx,z-l)(^x,; 



?cub 



?cub 



-icub 



+ \-)a 3 {-\® \f/>x,z+l){>l>x,z\], 

^? = £[|+m+|®|^-m)(^M 

x,z 

+ |-)<T 1 (-| ® 1^ + 1,^(^,2 |], 

x,y,z\ 

x,y,z 

+ h)a 3 (-| ® \^x,y,z+l)(^x,y 7 z\], 

^2i\^)c7 1 (M ® \^x-l,y,z)(^x,y,z\ 
x,y,z 

+ |-)cti(-| ® l^x,y+l, z )(^x,j/, z |], 



(12a) 
(12b) 
(12c) 
(12d) 
(12e) 



where position state \^ x ,y,z) = \^x) ® One 
complete step of the DTQW on a square and cubic lattice 
using two-state particle composes of the evolution in one 
axis followed by the evolution in the other, that is, 



W sq (9) = W* q (9)W* q 3 (9) 
W™ b (0)W™ b (0)W™ b (0), 



(13a) 



W cub (9) 
where W£(0) = SFJB^O) 

form of the evolution operators are, 



(13b) 

lz\ and W^" b (6») = 



z], respectively. Equivalent 



W^(^) = l x ® W CT3 (#) = ljf ® e-^3(fka) 



Wg(6) = W ai {9) ®l z = e~ iH ^ 6) ® l(J4b) 
Tv^ 3 " 6 (6>) = l x ® ly ® W<x 3 (0) = lx ® ly ® e- iff <-3(^4c) 
W™ b {6) = W ai (6) ® l y ® l z = e- iif 'i' 9 ' ® l y ® l(£4d) 

cub/ 



-iH a2 (6) 



The operator Wo- 3 (0) and ff^ (0) are given by Eq. Q and 
Eq. (l8|). Similarly, the operator for X and F axis, 



e - ip x cos(9) + sin(0) 
[cos(<9) + sin(<9) 

cos(0) - sin((9) 
- cos(#) + sin(<9) 

W ff2 (0) = S* 2 [B a2 (0) 

' cos(0) + sin(0) 
z cos(0) + i sin(0) 

cos(0) — sin(0) 
-zcos(0) + sin(0) 



g> 1] = e- iH "i^ 

cos(0) — sin(0) 
cos(0) — sin(0) 

- cos(0) - sin((9) 
cos(0) + sin(0) 

g> 1] = e - iH ^ e ) 

—i cos(0) + i sin(0) 
cos(0) — sin(0) 

z cos(0) + z sin(0) 
cos(0) + sin(0) 



(15a) 

(15b) 
(15c) 

,(15d) 



where Px and Py are the momentum operator on X and 
Y axis. The eigenvalues, A^ of W ai (0) and A^ of W (T2 (0) 
are same as the eigenvalues A^ with only a replacement 
of Px and Py in place of Pz in Eq. ([5]). Therefore, using 
the eigenvalues, eigenvectors and inverse of eigenvector, 
the Hamiltonian form for the evolution in X and Y axis 
are 



4 



In 



m 



2 v /cos 2 (6>) cos 2 (P_x) 

'"(f) 



2 A /cos 2 (6>) COS 2 (Py) 



cos(#) sin(Px) — i cos(Px) sin(0) sin(0) sin(Px) 

- sin(0) sin(Pv) cos(#) sin(P x ) + i cos(P x ) sin((9) 

cos(#) sin(Py) — i cos(Py) sin(0) —2 sin(0) sin(Py) 

-i sin(Py) sin(6>) cos(#) sm(P Y ) + z cos(Py) sin(0)_ 



• a\ (16a) 



• (J2 (16b) 



When = 0, that is, in absence of a coin operation A^ = 
e TiPz , A^ = e TiPx , and X Y = e TiPv . This reduces the 



Hamiltonian form to, H a3 (0) = 



Px 
. Px 



•cri, and H a2 (0) 



Pz 
P z 
P Y " 
Py 



-(73, ^(0) = 



•(72, respectively. 



If the initial state of the particle on a square lattice 

\c) = ^ 2 [\ i>+<it>],i* ln ) = ^niK>®i^>]' 

the state after t step [(W sg (<9)) £ ], 



and 



E E KUu> 



* (2) i 



t> ®|V>x,z>.(17) 



When = 0, t ) and t ) are given by the coupled 
iterative relations 



v( 2 ) 



*(x,y,t) 



l (x+l,y+l,t-l) 



\x+l,y-l,t-l) 



* (2) - a {2) 

*(x-l,y+l,t-l) "(a-l.y-l.t-l) 



(18a) 



(2) 

(x,y,t) 



*(x+l,y+l,t-l) 



a 



(i) 

(x+l,y-l,t-l) 



+a (2) + a (2) 

' "(a;-l,2/H-l,t-l) ' "(cc-l^-ljt-l) 



. (18b) 



In Fig. |2(a)| the probability distribution of the 50 step 
DTQW on a square lattice using the Pauli basis scheme 
without the coin operation (0 = 0) is shown. This prob- 
ability distribution obtained is identical to the ones re- 
ported for Grover walk on a four-state particle [33] and 
for the alternative walk on a two-state particle with ini- 
tial state, A=[\ \) + i\ t)] using Hadamard operator 
as the coin operation [36]. Each step of Grover walk 
is on a 2D is realized using the Grover diffusion oper- 
1111" 

I- 111 

II- 11 

III- 11 



ator G = - 



as coin operation Fol- 



lowed by S G = ^ 

| ® \^x-l,z+l)(^x,z\ + I 

)(-► I® I^E+M+lK^a;,*! 

state, \*g) = U\ I) -I 



|)(4|(8)|^-M-i>(^| + |t>(t 

<")(<- | ® + | 

on a particle in a specific initial 
> - | <-) + | -)►)]. The state after 




(a) 




(b) 



FIG. 2: Probability distribution of a 50 step two-state particle 
DTQW with the initial state |tf in ) = -^=[| t>] <S> |V>* > ® |^ ) 

on a square lattice using basis state of different Pauli operator for 
each axis (<xi for X axis and (73 for Z axis), (a) The distribution 
is after the evolution without the coin operation in both the axis 
(0 = 0) and same distribution is obtained for Grover walk, (b) The 
distribution is after the evolution with the coin operation 9 = tt/12 
in both the axis. 



t step of the Grover walk [S G (G <8> 1)*], 



i*?>= E E Kl t) u)+/5( ( l t) it>+C, t )i 

x= — t z=—t 

where /3(x, £)'s are given by the quadrupled iterative 



5 



relation coupling the X and Z axis 



/3< 4 > (*,*,*) 



B {1) + B {2) 

^(cc+l,z+l,t-l) ' ^(x+l,z+l,t-l) 



+s (3) + 5 (4) 

' ^(rc+l,^+l,t-l) ^ ^(rc+l,^+l,t-l) 



(20a) 



- 5 (2) 

^(X+1,2-I,t-1) ^(£C+l,Z-l,t-l) 



+ 5 (3) + 5 (4) 



(20b) 



8 {1) + 5 (2) 



-3 i3) + 5 (4) 

^(a:-l,«+l,t-l) ~ >^(x-l,z+l,t-l) 



(20c) 



+ B {2) 

^(x-l,z-l,t-l) ' ^(a-l.z-l,*-!) 



+s (3) - 5 (4) 

™(a;-l,2f-l,t-l) ^(aj-l,2;-l,t-l) 



(20d) 



From Eqs. (18) and Eqs. (20) we can note that for both, 



two-state walk and the Grover walk, the amplitude at 
any position (x, z) for a given time t is dependent on the 
amplitude at the four diagonally opposite sites at time 
t — 1. Therefore, starting from a specific initial state 
of two-state and four-state particle as discussed in this 
section, these amplitudes returns the same probability 
distribution. 

Unlike the Grover walk which is very specific to the ini- 
tial state and the coin operation, probability distribution 
with the two-state walk using the Pauli basis can be con- 
trolled by introducing the coin operation (0 ^ 0) and/or 
using different initial state of the particle. In Fig.|2(b) the 



probability distribution of the 50 step DTQW with coin 
operation, = n/12 is show to squeeze the distribution 
towards the diagonal of the square lattice. 

The basis states of as can be written as a superposi- 
tion of a basis states of the o\. Therefore, even in absence 
of a coin operation, [W^(0)WJ|(0)]* evolves the parti- 
cle in superposition of position space and implement a 
DTQW on a square lattice. Similarly, for a two-state 
walk on a cubic lattice using different basis states, due 
to the relationship between the basis states of the Pauli 
operators, the particle evolves in superposition of posi- 
tion space even in absence of coin operation. However, 
the coin operation can be effectively used to control the 
dynamics and the probability distribution of the walk. 



B. Triangular lattice 

The triangular lattice structures shown in Fig. [3] also 
has a three axis of propagation, X, Y, and Z. Therefore, 
the walk can be quantized using the eigenstates, 
and \—)a a of the Pauli operators a a where a = 1, 2, and 
3, as translational basis states. 

Unlike the cubic lattice the three axis for propagation 
in triangular lattice are not orthogonal to each other and 
hence, the evolution in one axis alters the evolution in 
other two axis as well. Therefore, the shift operator has 
to be defined according to the lattice structure. Shift 
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FIG. 3: Triangular lattice structure with labeling of lattice posi- 
tion in all the three axis, X, Y, and Z of propagation quantized by 
the basis states of the Pauli operators <ji, <J2 and (73, respectively. 




FIG. 4: Scheme for evolution of DTQW on: (a) The triangular 
lattice, starting from the middle, the arrow marks show the shift in 
position space during one step of DTQW evolution (evolution in Z 
axis followed by the evolution in the Y and X axis), (b) Kagome 
lattice structure with two axis of propagation at each lattice site. 
From lattice sites o, p, and q, we can see that they are associated 
with different combination of quantization axis. Starting from po- 
sition p, the arrow marks show the shift in position space during 
one step of DTQW evolution. The final positions are encircled. 



operator in each axis can be defined such that the unit 
shift in the main axis is accompanied by half of the unit 
shift in the other two axis. In Fig.[3j the labeling of the 
position is shown and for convenience we choose the two 
unit position shift in the main axis and one unit position 
shift in the other two axis. Therefore, the shift operation 
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for CT3 basis states is, 

S a3 = ^ U)U | ® \ip(x + l,y-l,z- 2)){iP(x, y, z)\ 

x,y,z 

+ J2 1 1) (t I ® \1>(x - 1, y + l,z + 2))(iP(x, y, z) [21) 



and the effective evolution operator 



iP 3 o\„(Q\„-iPa 

sin((9)e +iP3 cos((9)e +ip3 



cos((9)e- iP3 sin(0)e 



(22) 



where P3 = —Px ® Py ® 2Pz- Therefore, the Hamilto- 
nian H' a3 (0) will be in the same form of Eq. ([8| with a re- 
placement of Pz by P3. Similarly, the evolution operator 
W' ai {0) and (0), and the Hamiltonian form in the X 
and Y a xis, (0) a nd H' ao (0), will b e in t he same form 
of Eqs. ( |15b ), ( 15d| ) and Eqs. ( |16a ), (16b), respectively 
with the replacement of Px by Pi = 2Px ® Py ® —Pz 
and Py by P 2 = -Px ® ^Py ® JPz- Each Hamilto- 
nian, H' ai (0), H f a<2 {0), and H' a3 {0) evolve the state en- 
abling the interaction between the three quantization 
axis. Therefore, each step of DTQW can be realized by 
the evolution using one Pauli basis followed by the other 
W tri {0) = W^(0)W^ 2 (0)W^ 3 (0) as shown in Fig, [gaj 
We should note that the three effective Hamiltonian lor 
the operators, H' ai (0), H f a2 (0), and H' (0) together com- 
mute and complete Hamiltonian for each step of DTQW 
on triangular lattice can be written as 



H\0) = H'(9) + H'(6) + H'(8). 



(23) 



The choice of the order of the basis in which the particle 
is evolved is purely conventional for the triangular lattice. 
Even when = 0, due to the interplay between different 
Pauli basis for translation in each axis, a two-state parti- 
cle evolve in superposition of position space resulting in 
a diffused probability distribution. However, a coin op- 
eration with different for each axis can be extensively 
used for the evolution to get addition freedom to control 
the evolution a nd obt aine d the desired probability distri- 
bution. In Fig. 5(a) and |5(b)| we show the probability 
distribution of a 20 step DTQW without a coin opera- 
tion [H'(0) = H' ai (0) + #; 2 (0) + K 3 (0)] on a two-state 
particle initially in state | \) and | t), respecti vely. We 
can see t hat the probability distribution in Fig. 5(a) and 
Fig. |5(b)] are not symmetric distribution in po sition space 
but are symmetric to each other. In Fig. |5(c) we show the 
symmetric probability distribution obtained by introduc- 
ing a coin operation with = it/ A for only one operation 
during each step evolution [H' ai (0) + H' a2 (7r /4) + H' a3 (0)]. 



C. Kagome Lattice 

Kagome lattice structure can also be labeled the same 
way as the triangular lattice. The evolution operator and 
its Hamiltonian form in each basis (H f ai (0), H f a2 (0), and 
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80 
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FIG. 5: Probability distribution of 20 step two-state particle 
DTQW on triangular lattice, (a) The initial state is |^i n ) = | i 
) <S> IV^q) ® IV^o) ® I i^zo) and the walk is evolved without a coin op- 
eration, (b) The initial state is |^i n ) = | t) ® |V^ ) ® IVV)) ® l^o) 
and the walk is evolved without a coin operation, (c) The initial 
state is |^i n ) = I i) IV^o) ® \^y ) IV^o) anc ^ the walk is evolved 
with one coin operation [H' ai (0) + H' a2 (tt/A) + H' a (0)] to obtain a 
symmetric probability distribution in position space. 



H' a3 {0)) will be in the same form as presented for tri- 
angular lattice. But, unlike triangular lattice which has 
three quantization axis at each lattice site, kagome lat- 
tice shown in Fig. |4(b)| has only two quantization axis 
with four direction of propagation for the walk at each 
lattice site. The two quantization axis at each lattice 



site is not the same for all lattice sites. In Fig. 4(b) 
lattice sites o, p, and q have axis X and Z (<ji and 
03), X and Y (ai and 02), and Y and Z (a 2 and 03) 
as quantization axis, respectively. Therefore, to imple- 
ment each step of DTQW in kagome lattice certain sim- 
ple order for using evolution operators in different axis 
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has to be followed. For example, if the initial position 
is p (as marked in Fig. |4(b)|) , each step of DTQW can 
be realized by W kag {6) =W ai (0)W^ 3 (O)W^ 2 (O) and the 
effective Hamiltonian form for each step of DTQW will 
be H'{6) = H' ai {6) + H' a3 (0) + H' a2 (6). The main con- 
sideration in choosing the first axis is, not to pick the 
evolution in the axis which is nonexistent in that initial 
position. 



IV. CONCLUSION 

In this paper we presented a new scheme for the evo- 
lution of DTQW on 2D and 3D lattices using a two-state 
particle. Our scheme used different Pauli basis states as 
translational eigestate in different axis and showed that 
the coin operation is not a necessary requirement to im- 
plement a walk on 2D and 3D systems but can be used 
as an additional degree of freedom to control the dynam- 
ics. We also discussed the Hamiltonian form of evolu- 
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